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with  the  natural  boundary  condition  uxan  = 0 on  3 ft. 

Above  a is  a positive  definite  matrix  for  each 

x € ft,  b is  a vector,  and  n is  the  exterior 

normal.  Throughout  we  assume  that  the  functions 

2 _ 

aij'  ^i»  c are  class  C (ft)  and  the  bounded  do- 

2+a 

main  is  of  class  C . Then  the  principal  eigenvalue 
X*  exists  and  is  real  and  all  other  eigenvalues  X 

satisfy  Re  i > X*,  see  [ 1 6 ] • Here  a characteri- 
zation of  the  principal  eigenvalue  is  derived  moti- 
vated by  ideas  from  stochastic  control  theory.  This 
characterization  was  derived  earlier  by  Donsker- 
Varadhan  [ 6 ] using  different  methods.  The  stochastic 
control  theory  interpretation  suggests  a possible 
numerical  scheme  for  determining  the  principal  eigen- 
value and  eigenfunction  for  (1) , and  this  interpre- 
tation also  gives  a new  characterization  of  the 
smallest  eigenvalue  for  Schrodinger 1 s equation.  Fi- 
nally we  study  the  behavior  of  the  principal  eigen- 
value Xe  as  e + 0 when  a = cl  both  for  the 
Neumann  boundary  condition  = o and  the  mixed 
condition  + a u = 0.  The  above  results  are  only 
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outlined  here.  Some  of  these  results  have  appeared  / e. } 

in  (10],  [12].  I ,CA»o-< 
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Accompanying  Statement 


This  paper -4^s  a summary  of  an  invited  talk  presented  at 
the  International  .Conference  on  Stochastic  Analysis,  April  10-14, 
1978,  at  Northwestern  University.  This  paper  will  appear  in 
the  proceedings  of  £i^at  conference.  The  talk  consists  of  two 
main  parts.  First,  we  presents.our  recent  work  on  the  characteri- 
zation of  the  principal  eigenvalue.  We  have  been  able  to  derive 
a new  characterization  of  the  principal  eigenvalue  for  second  order 
linear  elliptic  partial  differential  equations,  not  necessarily 
self-adjoint,  with  both  natural  and  Dirichlet  boundary  conditions, 
and  also  give  a new  alternative  numerical  method  for  calculating 
both  the  principal  eigenvalue  and  corresponding  eigenvector  in  the 
case  of  natural  boundary  conditions.  The  principal  eigenvalue. 


if  appropriate  sign  changes  are  made,  determines  the  stability 
of  equilibrium  solutions  to  certain  second  order  nonlinear  partial 
differential  equations.  The  corresponding  eigenvector  enables  one 
to  determine  the  first  approximation  of  the  solution  of  the  non- 
linear equation  to  variations  of  the  initial  condition  from  the 
equilibrium  solution.  These  nonlinear  equations  are  important  in 
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the  applications.  For  these  reasons  it  is  important  to  have 
these  characterizations  of  the  principal  eigenvalue  and  eigen- 
vector. 

Our  method  converts  the  determination  of  the  eigenvalue 
and  eigenvector  to  determining  the  solution  of  a stationary 
stochastic  control  problem.  This  latter  problem  is  solved  and 
from  it  a numerical  scheme  arises  naturally.  This  method 
appears  to  have  applications  in  solving  other  problems. 

Secondly,  we  report  recent  progress  on  determining  the 
asymptotic  behavior  of  the  principal  eigenvalue  for  some 
singularly  perturbed  eigenvalue  problems  as  a small  nuisance 
parameter  tends  to  zero.  The  principal  eigenvalue  is  the 
optimal  value  for  a singularly  perturbed  stationary  stochastic 
control  problem.  We  are  thus  able  to  determine  the  asymptotic 
behavior  of  the  optimal  value  of  certain  stationary  stochastic 
control  problems . ^ 
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THE  PRINCIPAL  EIGENVALUE  FOR  LINEAR  SECOND 
ORDER  ELLIPTIC  EQUATIONS  WITH  NATURAL 
BOUNDARY  CONDITIONS 


Charles  J.  Holland 


I.  INTRODUCTION 


We  study  the  properties  of  the  principal  eigen- 
value X*  for  the  eigenvalue  equation 

-V- <uxa>  + uxb  + cu  = Xu  in  ft,  . (1) 

with  the  natural  boundary  condition  uxan  *=  0 on  3ft 
Above  a is  a positive  definite  matrix  for  each 
x € ff,  b is  a vector,  and  n is  the  exterior 
normal.  Throughout  we  assume  that  the  functions 


main  is  of  class  C . Then  the  principal  eigenvalue 
X*  exists  and  is  real  and  all  other  eigenvalues  X 
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satisfy  Re  X j>  X* , see  Ii6 1 • Here  a characteri- 
zation of  the  principal  eigenvalue  is  derived  moti- 
vated by  ideas  from  stochastic  control  theory.  This 
characterization  was  derived  earlier  by  Donsker- 
Varadhan  [ 6 ] using  different  methods.  The  stochastic 
control  theory  interpretation  suggests  a possible 
numerical  scheme  for  determining  the  principal  eigen- 
value and  eigenfunction  for  (1) , and  this  interpre- 
tation also  gives  a new  characterization  of  the 
smallest  eigenvalue  for  Schrodinger ' s equation.  Fi- 
nally we  study  the  behavior  of  the  principal  eigen- 
value Xe  as  e + 0 when  a = el  both  for  the 
Neumann  boundary  condition  = 0 and  the  mixed 
condition  + a2u  * 0.  The  above  results  are  only 

outlined  here.  Some  of  these  results  have  appeared 
in  [10],  [12]. 


II.  CHARACTERIZATION 

Define 

F (♦#V)  = *xa$*  + c*2  + bT(4a)-1b<j>2-|  div  b $2 

-|[(2a)“1b-V^]T(2a) [ (2a)“1b  - V*]*2.  (2) 
Theorem  1.  X*  satisfies 

X*  ■ min  max  [f  F($,V)dx  + [ i(b*n)  <|>2ds]  (3) 

♦€♦  v€c1(«)  Ja  ha4 

where  ♦ is  the  set  of  functions  $ of  class 
C2(0)  n C(ff)  with  *2  > 0 in  IT  and 
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f 2 -1 

1 4>  *=  1.  Note  that  whenever  (2a)  b is  a gradient, 

then  V in  (3)  is  chosen  so  that  Vx  = (2a)~*b.  In 
particular  when  b = 0 (3)  reduces  to  Rayleigh-Ritz. 

As  we  remarked  earlier  Donsker-Varadhan  [ 6 ] 
utilizing  their  powerful  methods  established  the 
analogue  of  (3)  for  the  Dirichlet  problem.  Their 
method  can  be  extended  to  cover  this  case.  We  es- 
tablished the  formula  (3)  using  stochastic  control 
theory  methods  in  [12]  and  outlined  how  to  establish 
the  result  also  for  the  Dirichlet  problem. 

Let  us  outline  our  derivation  of  (3).  If  X*, 
u*  > 0 represent  the  principal  eigenvalue  and 
normalized  eigenfunction,  then  (1)  can  be  rewritten 
aa 

-tr  au*  + u*(b+d)  + cu*  = X*u*  (4) 

for  an  appropriate  vector  function  d.  Now  u*  > 0 
in  IT.  Define  u*  = exp  (-ip),  then  ip  satisfies  the 
equation 

tr  a*xx  ” ♦x^x  “ + c = X*  (5) 

with  the  boundary  condition  <Pxan  * 0.  We  recognize 
that  (5)  is  the  Hamil ton- Jacobi  equation  for  a 
stationary  stochastic  control  problem  and  that  it  can 
be  rewritten  as 

tr  at  + min  [<>  v + (b+d+v) (4a) (b+d+v)  + c]  ■ Xt 
v€Rn  (6) 

A similar  stochastic  control  theory  interpretation 


The  stochastic  control  problem  is  the  following 
For  each  Lipschitz  function  w#  let  p be  the  in' 


variant  measure  associated  with  the  Ito  stochastic 


d£  *=  w(5)dt  + a(5)dz(t),  oa 


with  reflection  in  the  direction  -an 


With  each  control  w associate  the  cost 


[(b+d+w)(4a)  (b+d+w)+c]p  dx 


Then  X 


when  w = -b  - d + 2(u*)  (u*a) 


then  the  Min  C(w) 


b-d+2u  (u  a) 


Now  the  representation  (3)  is  obtained  by  dis 


equation  and  boundary  condition 


V.  [ (pxa)  - p (d+w) 1 * 0 in  ft 
(<Pxa)  - p(d+w)T)-n  = 0 on  3Q. 

A little  work  shows  that  every  w 6 it  can  be 
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V*(u2(-b  + 2(Wxa)T)  - 0 in  ft 

(-bT  + 2Wxa)*n  = 0 on  3ft.  (11) 

Now  W minimizes 

I (V)  = f ((2a)”1b-v3)T(2a) ((2a)“1b-v3)u2dx.  (12) 

x x 

Substituting  (10)  into  (8)  and  using  the  fact  that 

W minimizes  (12),  one  obtains  the  representation  (3). 

This  concludes  the  outline  of  the  proof. 

Corollary  1.  Min  Max  **  Max  Min  in  (3). 

♦ Cl  (II)  cl  (17)  * 

We  shall  outline  a proof  of  this  fact  below.  Note 
that  the  interchange  of  Min  and  Max  gives  the  repre- 
sentation 

X**=  Min  y (V) 

V€CX  (I?) 

where  y (V)  is  the  smallest  eigenvalue  to 

-tr  a$xx  + M(V)^>  = y<*  in  ft,  (14) 

M(V)  - c + bT(4a)_1b  div  b 

-y[  (2a) -1b-Vx) T (2a) [ (2a) _1b-v£] , 

with  the  boundary  condition  4>xan  + y(b*n)<J>  =0  on 

3ft.  It  is  not  necessary  that  y(b»n)  ^ 0 to 
guarantee  existence  of  the  smallest  eigenvalue.  Now 
the  fact  that  Min  Max  = Inf  Sup  « Sup  Inf  follows 
from  the  fact  that  the  functional  in  (3)  is  strictly 
'concave  in  v if  we  require  the  additional  assumption 
that  V(0)  * o.  Convexity  in  the  argument  $2  follows 
from  the  following  lemma  which  is  the  obvious  modifi- 
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cation  of  Lemma  2.2  in  Donsker-Varadhan  [ 6 ] . Note 
that  if  q = <J>2,  then  J (qxaq£)q-1dx  = j4<J>xa<J>£. 

Lemma  1.  Let  <2r  = {p:  p 6 C2  (8)  n C1  (IT)  # p > 0 in 
m,  Pxan  = 0 on  38).  If 

_ f t -1  f 70 <Pa) 

<1  € then  (q  aqMq  ax  = -inf  — qdx. 

Jn  x x P€^  p 

2 

Convexity  in  <{>  , concavity  in  V is  sufficient 
to  guarantee  that  Inf  Sup  = Sup  Inf,  see  Theorem  4.1 
in  Sion  [17],  The  fact  that  Sup  Inf  = Max  Min 
follows  from  the  fact  that  p (V)  exists  and  Remark 
2.3,  p.  175  in  Ekeland-Temam  [ 7 ] . 

III.  SCHRODINGER'S  EQUATION 

In  [10]  we  considered  the  Schrodinger  equation 

ih  (15) 

in  bounded  domains  under  the  cases  of  periodic,  zero 
Dirichlet,  and  zero  Neumann  boundary  conditions.  In 
(15)  H = -(h2/2m)Ax  + K,  h is  Planck's  constant 
divided  by  2n,  and  K >_  0 is  the  potential  function. 
Let  us  discuss  only  the  Neumann  case  here.  We  look 
for  solutions  of  (15)  of  the  form 

f(t,x)  = <p  (x)  •exp(-iAh“1t)  in  8 satisfying  = 0 
on  38.  Then  $ satisfies 

_ H*  - X*.  (16) 

Let  X*  denote  the  smallest  eigenvalue  of  (16) 
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and  the  corresponding  normalized  eigenfunction. 

We  now  formulate  a stochastic  control  problem  for 
1*#  $*•  Define  = exp(-Q/h) . Then  Q satisfies 
the  equation 

i . .*> 

(17) 


5^-  AxQ  + min(Qxw  + ^m|w|2  + K)  = 
w 


For  each  control  function  v let  pv  be  the  in- 
variant density  pv  generated  by  the  diffusion  pro- 
cess 

= v(£)dt  + (h/m)1/2dz(t)  (18) 

with  reflection  in  the  interior  normal  direction. 


Associate  with  each  control  v,  the  cost 

,2 


C (v) 


|at^|v| 


+ K]pv  dx. 


(19) 


Theorem  2.  X*  = MinvC(v)  where  the  minimum  is  taken 
over  Lipschitz  functions  v:  IT  -*•  Rn.  The  minimum  is 
obtained  when  v = (h/m)  (<J> */$*)  and  pv  = (<p*)2. 

Note  that  v plays  the  role  of  "velocity"  of  the 
process  (18)  . Thus  X*  can  be  thought  of  as  the 
minimum  energy  of  (19)  and  (<J>*)2  as  the  invariant 
probability  distribution  of  the  process.  This  inter- 
pretation corresponds  to  the  standard  physical  inter- 
pretations. See  Nelson  [15 J for  another  approach  to 
(15)  using  stochastic  differential  equations.  Of 
course,  the  interpretation  (18) , (19)  for  the  self- 
adjoint  problem  H*  - X*  is  Rayleigh-Ritz  in  disguise. 
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IV.  NUMERICAL  APPROXIMATIONS 


The  control  theory  interpretation  in  Section  2 
also  suggests  an  alternative  numerical  method  for 
approximating  the  principal  eigenvalue  and  eigen- 
function for  the  Neumann  problem  (1) . The  procedure 
is  to  discretize  the  verification  equation  (6)  ob- 
taining the  verification  equations  for  stationary 
control  of  a Markov  chain.  The  Markov  chain  problem 
is  the  discrete  analogue  of  the  stochastic  control 
problem  associated  with  (6) . The  Markov  chain  control 
problem  has  a special  form  which  can  be  solved  effi- 
ciently using  an  iterative  method  due  to  White  [19] . 

See  [10] , Section  3 for  details.  We  have  done  some 
numerical  work  on  simple  examples  using  this  method, 
however,  the  convergence  of  this  method  has  not  been 
justified  theoretically. 

This  type  of  numerical  approximation  can  also  be 
used  for  the  periodic  eigenvalue  problem.  However, 
the  formal  approximation  suggested  in  [10]  for  the 
Dirichlet  problem  fails. 

V.  SINGULARLY  PERTURBED  PROBLEMS 

Here  the  asymptotic  behavior  as  e + 0 of  the 
principal  eigenvalue  to  the  singularly  perturbed  eigen- 
value problem 


-eA$  - $x  (b  (x)  +eb  (x) ) + ( (c  (x)  +ec  (x) ) <J>  = A4>  (20) 
in  a bounded  convex  domain  ft  with  boundary  data 
+ a(x)4>  = 0 on  is  studied.  The  behavior  of 

A has  been  extensively  studied  for  the  Dirichlet 

problem  with  c = 0;  see  the  papers  [3]#  [9]/  [18] . 

There  the  behavior  can  vary  from  Ae  -*■  0 to  Ae  -►  «® 

as  e + 0 depending  upon  the  nature  of  the  vector 

field  b.  For  the  zero  Neumann  problem  such  a wide 

range  of  behavior  is  not  possible  since  one  knows  that 

min (c (x) +ec (x) ) Ae  <_  max  (c  (x)  +ec  (x) ) . However,  the 
x€?7  x€ft 

behavior  of  the  principal  eigenvalue  for  Neumann  and 
mixed  problems  is  of  interest,  especially  in  the 
stability  analysis  of  equilibrium  states  of  nonlinear 
diffusion  processes  occurring  in  chemical  reactors. 

See  H.  Keller  [14]  who  treats  the  eigenvalue  problem 
in  one  spatial  dimension  by  transforming  the  problem 
to  the  standard  self-adjoint  Sturmian  form. 

The  problem  with  the  zero  Neumann  condition 

du 

■jjj  *=  0 is  treated  first.  The  mixed  problem  with 
o | 0 will  then  be  solved  by  transforming  the  mixed 
problem  into  a Neumann  problem.  The  properties  of  the 
principal  eigenvalue  Ae  depend  upon  the  nature  of  the 
vector  field  b + eb,  as  is  evident  from  (7)*- (8).  We 
need  the  following  definition.  Throughout  let  £°(t) 
be  the  solution  to  the  differential  equation 


<U°(t)  = b(£°(t))dt,  e°(0)  = x € 17  (21) 

with  reflection.  Reflection  means  that  if  y € 3ft 

and  b(y)*n(y)  >_  0,  n the  exterior  normal  as  before, 

then  £y(t)  travels  along  the  boundary  with  velocity 

equal  to  the  projection  of  b on  the  tangent  space  to 

fl  at  £°(t).  Otherwise  £°(t)  is  the  solution  to 

(21) . This  prevents  the  path  from  escaping  from  ft. 

See  Bensoussan-Lions  [ 2 ) or  Holland  [11]  for  a more 

detailed  discussion  of  the  path. 

First,  the  self-adjoint  case  b = b is  disposed  of. 

Lemma  2.  If  b = b = 0,  then  Xe  -*•  min  c(x)  as 

x€ft 

c + 0. 

Proof:  This  follows  immediately  from  the  Rayleigh-  • 
Ritz  representation. 

Theorem  3.  Suppose  that  the  equation  (21)  with  re- 
flection has  a globally  (for  all  x in  I?)  asympto- 
tically stable  rest  point  x.  Then  Xe  c(x)  as 
e 4 0. 

Remark . A rest  point  x is  a point  in  17  such  that 
5j(t)  « x for  all  t.  If  x € ft,  then  b(x)  * 0. 

If  x 6 !7,  then  either  b(x)  « 0 or  else  b(x)  is 
parallel  to  xv(x) . 

Outline  of  proof:  Let  Xc  denote  the  principal  eigen- 
£ 

value  and  $ >0  the  corresponding  eigenfunction. 

Define 
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uE(x,t)  *=  4>c  (x) e”*Et. 


(22) 


Then  u satisfies 

eAu  + u (b+eb)  - (c+ec)u  - u 


x' ~t 

dll 

with  the  boundary  condition  = 0 on  aftxIO,®). 
Then  using  the  Ito  differential  rule  we  have  that 


(23). 


ue(x,t)  = E*(CE(t))expft-c(d(s))ds 

X J n X 


(24) 


where  £^(t)  is  t^ie  s°luti°n  to  tile  Ito  equation 


<UE(t)  = (b+eb)  (CE(t))dt  + (2e)1/2dz(t)  ,£E(0)  = x. 


(25) 


z(t)  Brownian  motion  as  before,  with  reflection  in  the 
interior  normal  direction  at  the  boundary  of  See 

I 2l  for  a construction  of  the  process. 

The  proof  is  then  completed  using  the  represen- 
tations (22)  and  (24)  for  sufficiently  large  t.  This 
type  of  approach  follows  ideas  of  Kac  [13].  The  proof 
requires  the  estimate  that  in  convex  domains 


Pr  [||5*C)-e°(-)|l  > p]  i 0 
0<t<t ' x x 


as  c 4 0 for  any  fixed  y > 0.  This  estimate  is  de- 
rived in  [11]  (see  equation  (15)  there)  and  depends 
upon  estimates  (3.13),  (3.14)  in  [ 2]  for  reflecting 
processes  in  convex  domains. 

Example . If  ft  is  the  unit  disc  in  R2  and 
b - (bj^bj)  - (1,-x2) , then  X€  •*>  c(l,0)  as  c * 0. 
Theorem  4.  Suppose  that  the  differential  equation  (21) 
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I 


with  reflection  has  a limit  cycle  r whose  domain  of 

attraction  consists  of  all  of  Ti  except  for  isolated 

critical  points  in  II.  Further,  at  a critical  point 

x*  assume  that  the  eigenvalues  of  bx(x*)  have  only 

positive  real  parts.  Then 

rT* 

j c(£°(t)  )dt  as  e + 0 

where  T*  is  the  period  of  r and  x is  an  arbi- 
trary point  on  r. 

Proof:  The  proof  is  similar  to  the  proof  of  Theorem  3 
except  one  has  the  added  complication  of  possible 
critical  points  x*. 

We  now  consider  the  mixed  problem  for  the  cases 
considered  in  Theorems  3 and  4. 

Theorem  5.  Let  H be  a smooth  solution  to  AH  - H * 0 
in  ft  with  — a on  3ft.  Suppose  b satisfies 
the  assumptions  of  Theorem  3.  Then 

Ac  ♦ c(x)  + VH (x) *b (x)  as  e + 0.  Suppose  b satis- 
fies the  assumptions  of  Theorem  4.  Then 

fT* 

* fr  j c (C-(t) ) + VHUf(t)) -bu£(t))dt  as  e + 0. 

Proof:  The  result  follows  by  considering  the  eigen- 
value problem  for  v - eH^u. 

Example . The  limits  for  Ac  in  Theorem  5 earn  be 
evaluated  in  the  following  case.  Suppose  8 c Rn  and 
b satisfies  the  assumptions  of  Theorem  3 with 

I 


i 
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x € 3ft.  Then  b(x)  = |b(x)|n(x)  and  hence 
Xc  c(x)  + <x(x)  jb(x)  | . Thus  if  ft  = (0,1) , b > 0 

on  [0,1],  then  Xe  -►  c(l)  + a(l)b(l)  as  e + 0. 
This  last  remark  is  the  content  of  the  limit  result 
(1.7)  in  [14]. 


ACKNOWLEDGMENT 


The  author  would  like  to  thank  S.  R.  S.  Varadhan 
for  some  stimulating  discussions. 


REFERENCES 


1.  Anderson,  R. , Orey,  S.,  Small  random  perturba- 
tions of  dynamical  systems  with  reflecting 
boundary,  Proc.  Symposium  Pure  Math.  Vol.  XXXI, 
Univ.  111.,  1976,  pp.  1-4,  American  Math  Soc. 

2.  Bensoussan,  A.,  Lions,  J.,  Diffusion  processes 
in  bounded  domains  and  singular  perturbation 
problems  for  variational  inequalities  with 
Neumann  boundary  conditions  in  Probabilistic 
Methods  in  Differential  Equations,  Springer 
Verlag  Lecture  Notes  in  Mathematics  #451(1975), 
8-25. 

3.  Devinatz,  A.,  Ellis,  R. , Friedman,  A.,  The 
asymptotic  behavior  of  the  first  real  eigenvalue 
of  second  order  elliptic  operators  with  a small 
parameter  in  the  highest  derivatives  II,  Indiana 
Univ.  Math.  J.  23(1974)  991-1011. 

4.  Devinatz,  A.,  Friedman,  A.,  Asymptotic  behavior 
of  the  principal  eigenfunction  for  a singularly 
perturbed  Dirichlet  problem,  preprint. 

5.  Donsker,  M. , Varadhan,  S.  R.  S.,  Asymptotic 
evaluation  of  certain  Wiener  integrals  for  large 
time,  in  Proceedings  of  the  International  Con- 
ference on  Integration  in  function  spaces. 
Clarendon  Press,  Oxford  (1974)  82-88. 

€.  Donsker,  M. , Varadhan,  S.  R.  S. , On  the  principal 
eigenvalue  of  second  order  differential  operators. 
Comm.  Pure.  Appl.  Math.  29(1976)  595-621. 


14 


I 


Ekeland,  I.,  Teman,  R. , Convex  Analysis  and 
Variational  Problems.  North-Holland  (197  6)  . 
Fleming,  W. , Exit  probabilities  and  optimal 
stochastic  control,  preprint. 

Friedman,  A.,  The  asymptotic  behavior  of  the 
first  real  eigenvalue  of  a second  order  elliptic  • 
operator  with  a small  parameter  in  the  highest 
derivatives,  Indiana  Univ.  Math.  J.  22(1973) 
1005-1015. 

Holland,  C. , A new  energy  characterization  of  the 
smallest  eigenvalue  of  the  Schrodinger  equation 
Comm.  Pure  Appl.  Math.  30(1977)  755-765. 

Holland,  C.,  The  regular  expansion  in  the 
Neumann  problem  for  elliptic  equations  Comm,  in 
Partial  Differential  Equations  1(1976)  191-213. 
Holland,  C. , A minimum  principle  for  the  smallest 
eigenvalue  for  second  order  linear  elliptic 
equations  with  natural  boundary  conditions.  Comm. 
Pure  Appl.  Math.,  to  appear. 

Kac,  M.,  On  some  connections  between  probability 
theory  and  differential  and  integral  equations, 
Proc.  of  the  2nd  Bekeley  Symposium  (1950) 
pp.  189-215. 

Keller,  H.,  Stability  theory  for  multiple 
equilibrium  states  of  a nonlinear  diffusion 
process,  SIAM  J.  Math.  Anal.  4(1973)  134-140. 
Nelson,  E. , Dynamical  Theories  of  Brownian 
Motion , Princeton  University  Press,  1967. 

Protter,  M. , Weinberger,  H. , On  the  spectrum  of 
general  second  order  operators.  Bull.  Amer.  Math. 
Soc.  72(1966)  251-255. 

Sion,  M.,  On  general  minimax  theorem,  Pac.  J. 
Math.  8(1958)  171-176. 

Ventsel,  A. ,Freidlin,  M. , On  small  random 
perturbations  of  dynamical  systems,  Russ.  Math. 
Surveys  25(1970)  1-56,  Uspekhi  Math.  Nauk, 
28(1970)  3-55. 

White, D.,  Dynamic  programming,  Markov  chains, 
and  the  method  of  successive  approximations, 

J.  Math.  Anal.  Appl.  6(1963)  373-376. 


